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We show non-perturbatively that the power spectrum of a self-interacting scalar field in de Sitter 
space-time is strongly suppressed on large scales. The cut-off scale depends on the strength of the 
self-coupling, the number of e-folds of quasi-de Sitter evolution, and its expansion rate. As a conse- 
quence, the two-point correlation function of field fluctuations is free from infra-red divergencies. 
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A central building block of our current understanding 
of the Universe is the idea of cosmological inflation. It can 
be realized by a single scalar field, called inflaton. The 
quantum fluctuations of this field and of the space-time 
seed today's structures in the Universe. This is well un- 
derstood in the case of small fluctuations. However, some 
models of inflation imply the existence of large quantum 
fluctuations, e.g., the so-called chaotic scenario [T], which 
also leads to the idea of eternal inflation [2]. 

Our task is to study the implications of large quan- 
tum fluctuations on structure formation and the evolu- 
tion of the Universe. Of special importance is the power 
spectrum, which is observable via the cosmic microwave 
background radiation [3]. 

A toy model to understand the physics of large quan- 
tum fluctuations is a self-interacting scalar field $ in de 
Sitter space-time. In the traditional approach, based on 
perturbative quantum field theory in curved space-time, 
the multi-point correlation functions of $ generically ex- 
hibit infra-red divergencies (c.f. 

Related to these issues is the presence of large fluctua- 
tions over super-horizon distances, which can invalidate 
the entire semi-classical approximation [S]. This break- 
down of semic-classical methods is due to the fact that 
higher loops are not suppressed, meaning that all orders 
in the loop expansion have comparable sizes. Hence, non- 
peturbative methods are needed for a full description and 
deeper understanding of inflationary cosmology. 

Stochastic inflation [6Hl4] provides one of the very few 
of such approaches. Its idea lies in splitting the quantum 
fields into long- and short-wavelength modes, and view- 
ing the former as classical objects evolving stochastically 
in an environment provided by quantum fluctuations of 
shorter wavelengths. Given the de Sitter horizon c/H 
as a natural length scale of the problem, one then fo- 
cusses on the "relevant" degrees of freedom (the long- 
wavelength modes) and regards the short-wavelength 
modes as "irrelevant" ones, where "short" and "long" 
are subject to the horizon. 

The most simple setup provides a fixed cosmological 
background, in which the dynamics of a scalar test field 
$ is analyzed. If $ is free, massive and minimally cou- 
pled, one obtains after splitting into long and short wave- 
lengths, $ = -|- 0, an effective equation of motion of 



generalized Langevin-type, 

(□ + M^) <^ h. 



(1) 



The quantity h is a Gaussian-distributed random force 
with zero mean. 

In [131 [H] we applied replica field theory together with 
a Gaussian variational method to stochastic inflation. 
Extending early studies, that mainly focussed on homo- 
geneous flelds and thus restricting attention to the time 
evolution of $, we presented a method to calculate arbi- 
trary two-point correlation functions. 

In this work we extend our previous results [T31 [H] 
to include self-interactions of a scalar fleld in de Sitter 
space-time. For the speciflc example of a quartic self- 
interaction we calculate the power spectrum and show 
that self-couplings cause a damping of this quantity on 
large scales. This therefore solves the problem of infra- 
red divergencies of two-point correlation functions. Fur- 
thermore, our results are independent of the precise way 
of splitting into long- and short-wavelength modes and 
hence do not suffer from any related ambiguity. 

As a starting point, we use the Lagrangian 



C 



(2) 



where 4> is a massless, minimally-coupled, real scalar field 
with quartic self-coupling constant A. Greek indices run 
from to 3. We assume a de Sitter background geome- 
try, {g^u) = diag(l, —a{ty,—a{t)'^, —a{t)^) with the scale 
factor a{t) = exp{Ht). For convenience we use h = c = 1. 

Let $Q be a free field, being subject to ^ with A = 0. 
It might be decomposed as 



$o(i, fc) = a(fe) Uoit, k) + h.c, 



(3) 



with the modulus of the comoving momentum k :— \k\. 
The annihilation and creation operators, a and a^, obey 
the usual commutation relations. 

In terms of conformal time r, with a(T) = (iJr)^^, the 
rescaled mode functions Vq{t, k) :— a(r) Uq{t, k) fulfil the 
mode equation 



= 0, 



(4) 
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where primes denote derivatives with respect to r. So- 
lutions to Q are fixed by requiring that for very short 
wavelengths the effect of space-time curvature becomes 
irrelevant, and thus a plane-wave solution should be 



where h is a Gaussian-distributed random variable with 



obtained, limfc/^-foo ■'^o(''"' ^) — e /v2fc. The factor 
l/v'2 k is fixed by the canonical commutation relations 
of $0 ^i^d its conjugate momentum. At late times, the 
leading term of the solution to Q reads 



zto(A;«l/|r|)^- 



iH 
V2W 



(5) 



An object of central interest in cosmology is the power 
spectrum V{k). Its relation to the propagator 

G{k){2Tr'fS^{k-k') = $(fc')|f^), (6) 

where |r2) is the Bunch-Davies vacuum, is given by 

fc3 



V{k) 
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2 Gik). 



(7) 



On supcrhorizon scales (A; <C 1/|t|) one finds for the free 
massless case (subscript "0" ) a scale- invariant spectrum: 
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To go beyond Equation ([8|, we first split the field $ 
into a short- and a long-wavelength part, $ = -f 1^9, 
where we use the filter function F^;, specified by its deriva- 
tive: 



Kiy) = 



'0 

iV exp ( 1 - 
.0 



1 - 



V < -K, 
y e [-k,k], 
y > +K, 



(9) 



with y := fc|r| — e, cutting out wave numbers below e/jrl 
with a cutting width k. In the limit k — > 0, F^; ap- 
proaches the step function 8(fc|r| — e)- The constant 
N := c^tT V2/(5.3k) in ^ is a normalization factor, 
ensuring J dyF^{y) = 1. Throughout this work we choose 
K = 10~^ and e = 10^^, although our main statements 
are virtually independent of these quantities. 

Having introduced the precise way of splitting into 
long- and short-wavelength modes with the filter func- 
tion (|9]), we now consider the form of the induced noise 
terms. For general self- interactions, they are non-linear 
in the short-wavelength modes. However, if one is inter- 
ested in the late-time behavior, or more precisely in the 
leading-ln(a(i)) contribution, one may restrict to linear, 
Gaussian-distributed noise terms. This has been argued 
already a long time ago by Starobinsky [5] and has been 
rigorously proven by Woodard 

The stochastic field equation for model ^ reads then 



Dip + AXif^ =h, 



(10) 



h = 0, h2 



(11) 



where A is a known function, depending on derivatives of 
the mode functions in ([s]), see e.g., [5] for a more detailed 
presentation. The "bar" in (11) and below denotes the 



average over the noise due to quantum fluctuations of 
short wavelengths. 

Different filter functions have been intensively dis- 
cussed in [13]. In [13] we showed for free fields that filter 
functions with compact support allow us to avoid infra- 
red divergencies. Below we show that they are also absent 
if quartic self- interactions are turned on. In fact, this re- 
sult in independent of the choice of a filter function. 

Let us now briefly summarize the program we will per- 
form next: As we described in detail in [131 [M], we first 
Wick-rotate to Euclidean signature and use the replica 
trick jl5] . Then we introduce a suitable variational action 
and determine its form (especially its replica structure) 
from a Feynman- Jensen variational principle |16| . This 
will allow us to go beyond ordinary perturbation theory 
(c.f. [13l [TU [17] ) and to obtain an analytic expression for 
the full power spectrum. 

After Wick-rotating we proceed with the replica trick. 



6}{xi) . ..Si{Xn) 



1 (5" 

= hm ——7 — r — 

m-yO m 0](Xi) . . 



ln(2"[j]), 



(12) 



where Z[j] is the generating functional depending on an 
external current j. m denotes the number of replicas, 
labelled by the indices a,b, . . . . Furthermore we define 
the replicated action 5^™^ via 



/m / m \ 

n^^bJexp -^5[^„j] 

a=l \ b=l / 



(13) 



n2?bJexp(-5M[{^},j] 



Besides terms diagonal in replica space (cx Sab), it also 
contains the non-diagonal part 

(14) 

originating from the average over noise. 

We apply the Feynman-Jensen variation principle and 
therefore define a Gaussian variational action 

^irHM]:^;^^ / ^a{t,k)G-\b{t,k)Mt,~k), 

(15) 
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with := Jdt and := Jd^k/{2'K)^. We make the 
ansatz for the inverse propagator 



Gn 



6n 



(16) 



The self-energy matrix \(j5ab — (^ab\ mimics the diagonal 
and the non-diagonal parts in (13 1, respectively. 



Maximizing the right-hand side of the Feynman- Jensen 
inequality 

ln(Z)>ln(Z_) + ^5i7)-5(™)^ , (17) 
wherein the subscript "var" refers to the variational ac- 



tion ( 14 ) , yields the replica symmetric solution 



771 
AJab, 



Tr[G,J, 



(18a) 
(18b) 



with a ^ h. The 777 x r77-matrix J is defined by JJafc = 1 
for all a, b. 



To solve the implicit equations (18a) and (18b) we in- 



vert ( 16 ) by means of an expansion in the number of 



replicas 771. At leading order we find 

-1 



Gab 
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(19) 



From this quantity we extract the full physical propaga- 
tor G{t,k) via (c.f. HTj) 



G{t,k) - lim -Tr|G„b(t,fc)l 

m-!-0 777 



(20) 



Its explicit form is rather lengthy and shall not be given 
here, but the corresponding power spectrum is plotted in 
Figure [1] At late times one can identify two regimes: 



r{t,k) 




: k < fc^(t), 

: fc,(i)<fc<l/|T(<)|. 

(21) 



The comoving wave number k^ at which the large-scale 
behavior of 7^(i, k) changes significantly is determined by 
tj from ( 18a[ ). As a function of the number N of quasi-de 

Sitter e-folds, for TV at least a few and with N{t = 0) = 0, 
k^ reads 



27r2 



(22) 



The factor in front of the curly brace in equation ( |2l[ ) is 
the standard value of the scale-invariant power spectrum. 
The large-scale behavior \k <^ fc*(i)] ofV{t,k) follows 
from the fact that: a) the quantity a is /c- independent 
[c.f. (18a)], b) the second term in (191 is subdoniinant 



compared to the first [hence Gab{t, fc — )■ 0) = const.], and 
c) the relation of 7-'(t, k) to G(t, k) involves a factor k^ . 
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Figure 1: Power spectrum 'P(fc) of a massless test field with 
quartic self-coupling as a function of comoving momentum. 
Upper panel: Results for various values of A, with A'^ = 6. 
Lower panel: The same as above but for various numbers of 
e-folds, where the coupling has been fixed to A = lO"^'^. Also 
indicated is the value at which the power spectrum has 
decreased to half of its amplitude (normalized to 1 here). 



We observe in Figure [T] that the power spectrum is 
heavily suppressed on large scales in agreement with 
This damping becomes more pronounced as the 



(21) 



self-coupling A is increased and for a large number of e- 



folds [c.f. (22)]. Hence, the self-coupling breaks the scale 
invariance of V{t, k). 

On subhorizon scales one finds V{k ^ I/ItI) ~ k^ . 
This might be understood from: a) the fact that quan- 
tum effects in stochastic inflation only significantly mod- 
ify large scales, b) the behavior of the free mode function 
Mo(fc > l/|r|) (,) the factor k^ in 

The derived large-scale suppression solves the problem 
of infra-red divergencies of real-space correlation func- 
tions: While in a scale-invariant theory the two-point 
function diverges in real space, 



G{t,x) 



(27r)3 



'■G{t,k) cx 



dfc sin(fc) 
~k V~ 



— ?> 00, 
(23) 



the theory with correctly-resummed quantum effects will 
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be finite. For a complete understanding of quantum ef- 
fects in inflationary cosmology one would need to include 
metric fluctuations. 

Let us now study if a cut-off at could be observable. 
We assume quasi-de Sitter inflation and a sudden reheat- 
ing to the radiation-dominated Universe after N e-folds. 
This gives for the physical damping scale today, 



3 
2^ 



H 



Ho. (24) 



Hq is the present value of the Hubble rate, and T,„i, and 
To denote the temperatures at reheating and today, re- 
spectively. With A ~ 10-^3 and T„„ ~ 10^^ g^d 
for iV « 51 that today the physical damping scale is of 
the order of the horizon: 



Hn 



(25) 



For a much larger number of e-folds this cut-off is unob- 
servable. 

Other scenarios jTSl Hn] with a finite number of e-folds 
also lead to a cut-off in V{k). However here it is mainly 
the self-interaction, which is responsible for the large- 
scale damping. This can be easily seen by considering 
that the A = result is in our case ([s]). 

Suppression of the power spectrum in the infra-red 
could also influence the cosmic microwave background 
radiation. This issue was brought into focus by re- 
cent observations [31 |20l [21], which suggest a lack of 
power on the largest observable scales. Despite a cut- 
off in the primordial power-spectrum about the Hubble 
scale, the integrated Sachs- Wolfe effect [22] can regener- 
ate power on the largest observable scales in the cosmic 
microwave background. Mortonson and Hu |23j recently 
provided new upper bounds on such a cut-off. They 
found fc,„. < 5.2x 10-'^Mpc"^(95% C.L.) using polariza- 
tion data. This value is close to and well consistent with 
the above estimate of fcj'' ~ Hq ~ 2.4 x lO^^'Mpc^^ 

So far we considered tiny self-couplings. However, 
thanks to the non-perturbative nature of the used meth- 
ods, our result also applies for stronger couplings. One 
important example is that of the Higgs boson, which 
has A = 0(0.1) and — in the perturbative approach — 
also suffers from infra-red divergences in quasi-de Sitter 
spaces. In that case is much closer to the horizon scale 
during inflation. 

The occurrence of infra-red divergencies in standard 
perturbation theory is linked to a breakdown of semi- 
classical methods of quantum field theory in de Sitter 
space, which has been recently pointed in reference [S]. 
This is shown therein for a scalar field with quartic self- 
interaction. It has been found that this break-down can- 
not be cured by including any finite number of loops. 
It rather concerns the failure of the entire semi-classical 
calculation. 



To summarize, quantum effects in inflationary cosmol- 
ogy significantly modify the large-scale evolution of quan- 
tum flelds. Within the framework of stochastic inflation 
and using the methods of replica field theory, we have 
shown for the specific example of a self-interacting scalar 
field in de Sitter space-time, that the power spectrum is 
free from infra-red divergencies due to a large-scale cut- 
off. Our results are free from any ambiguity associated 
with the choice of a particular filter function. Further- 
more, our findings have been obtained non-perturbatively 
and are hence not plagued by a break-down of standard 
perturbation theory in de Sitter space-time (c.f. [5 ). 
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